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£N) ■ Halos are biased tracers of the dark matter distribution. It is often assumed that the initial 

patches from which halos formed are locally biased with respect to the initial fluctuation field, 
<^^> meaning that the halo-patch fluctuation field can be written as a Taylor series in the dark matter 

O^l ■ density fluctuation field. If quantities other than the local density influence halo formation, then 

this Lagrangian bias will generically be nonlocal; the Taylor series must be performed with respect 
to these other variables as well. We illustrate the effect with Monte-Carlo simulations of a model 
in which halo formation depends on the local shear (the quadrupole of perturbation theory), and 
provide an analytic model which provides a good description of our results. Our model, which 
C^*) \ extends the excursion set approach to walks in more than one dimension, works both when steps 

in the walk are uncorrelated, as well as when there are correlations between steps. For walks with 
correlated steps, our model includes two distinct types of nonlocality: one is due to the fact that 
the initial density profile around a patch which is destined to form a halo must fall sufficiently 
. steeply around it - this introduces fc-dependence to even the linear bias factor, but otherwise only 

affects the monopole of the clustering signal. The other is due to the surrounding shear field; this 
affects the quadratic and higher order bias factors, and introduces an angular dependence to the 
clustering signal. In both cases, our analysis shows that these nonlocal Lagrangian bias terms can 
Q ■ be significant, particularly for massive halos; they must be accounted for in analyses of higher order 

?— i 1 clustering such as the halo bispectrum in Lagrangian or Eulcrian space. Although we illustrate these 

ryj , effects using halos, our analysis and conclusions also apply to the other constituents of the cosmic 

■ web - filaments, sheets and voids. 
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I. INTRODUCTION 



' The virializcd halos which are identified in simulations of gravitational clustering are biased tracers of the underlying 



matter field. Typically, this bias is described in two ways, cither by relating the halo and mass fields at the time the 
halos were identified (e.g., the present), or by identifying the patches in the initial conditions which are destined for 
fN) . form halos, and describing the bias between these patches and the initial mass fluctuation field [l[. These are known 
as Eulerian and Lagrangian bias, respectively. In either case, the simplest models assume that this bias is local, 
meaning that the biased field can be written as a (deterministic) function of the mass field. However, the nonlinearly 
evolved mass field is a nonlocal function of the initial one, so Lagrangian and Eulerian bias cannot both be local |3|-|5|. 
' It has recently been noted that neither of the two best studied models of Lagrangian bias, peaks theory Q and the 
excursion set approach 0] , are local. This is because both approaches predict that the abundance of biased tracers 
(peaks or halos) should depend, not just on the local values of the overdensity field, but on derivatives of the field as 
well [H . This gives rise to a rather specific form for nonlocal bias, in which the bias is most naturally described in 
Fourier space, where it is fc-dependent, even at the linear level [9l4ll|. The main goal of the present work is to explore 
models in which the nonlocality of bias is qualitatively different, arising only at second order, and associated with 
anisotropies in the initial field. This source of nonlocality is generic to models in which halos form from an anisotropic 
collapse [lH , for which there is considerable evidence [l3| . 

Section |TT] describes the relation between Lagrangian and Eulerian bias in nonlocal models, and Section InTl describes 
our excursion set based treatment of the origin of nonlocal Lagrangian bias, providing an analytic description of the 
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effect. The technicaf problem which is the subject of this section is to provide an accurate formula for the first crossing 
distribution of a barrier by n-dimensional walks with correlated steps; we use Monte-Carlo simulations of such walks 
to illustrate the accuracy of our analytic formulae. A final section summarizes our findings. 



II. NONLINEAR, NONLOCAL HALO BIAS AND ITS EVOLUTION 

The excursion set model of halo abundances and evolution assumes that it is possible to identify those patches in 
the initial fluctuation field which are destined to form halos by a later time (e.g., the present) 0. In what follows, 
we discuss why, in this case, halo bias is generally expected to be nonlocal both in the initial conditions and at later 
times. 



A. The spherical evolution model: Local Lagrangian bias and evolution 

The simplest implementations of this approach assume that halos form from a spherical collapse. In this case, the 
initial overdensity of a patch plays an important role in determining whether or not it will form a halo (T3 |. Almost 
all studies which incorporate the spherical collapse model into the excursion set approach assume that the initial 
overdensity is the only parameter which determines halo formation (see Q for a recent exception). As a result, these 
studies find that halos are locally biased versions of the initial density fluctuation field If the environment of a 
halo is spherically symmetric, then it is natural to expect the spherical model to describe both the evolution of the 
patch which becomes a halo, and the evolution of its surrounding environment. Since the spherical evolution model 
yields a local deterministic mapping between the initial and evolved densities, the fact that halo bias was local with 
respect to the initial density field means that it remains local with respect to the evolved field. The bias factors are 
different of course, but they are easily calculated Q. 

Schematically, if an initial volume Vq is overdense by Sq, then, in the excursion set approach, the halo overdensity 
(averaged over all such spheres) is 



(Nm\S ) = v b\8§ 
n(m)V £t k\ 



1 + h(m\5a) = ^£^F = 1 + 2^ (1) 

fc>0 



where the numerator is an average over all cells of volume Vq and density So, and n(m) in the denominator denotes the 
number density of halos of mass m (i.e., the integral of the numerator over all allowed values of <5o)- The right-hand 
side is explicitly a function of Jo only, so the Lagrangian bias is local. The Lagrangian bias factors associated with 
the spherical model satisfy 

$c b k = v k ~ X %iW where v = 8 c /a(m), (2) 

and the H n are (the probabilist's) Hermite polynomials, and 8 C is the critical density required for spherical collapse. 
The Eulerian halo overdensity is defined similarly: 

i^&< \m .a - (N m \5 (M,V)) , M ^b k 8 k 

1 + 8 h (m\M, V) = = i + ^ (3 ) 

V ; k>0 

The Eulerian halo bias factors are obtained upon noting that M/pV = 1 + 5 and that the spherical model yields a 
local monotonic relation between 1 + 8 and <5oi so that the expression above can be written as a series in 8. Specifically, 
the Eulerian bias factors in the spherical collapse model are: 

bx = l+ V -^- = l + b\ (4) 

01 

b 2 = ^b\ + b\ (5) 
796 , T 13 lT ,, 

b 3 - -— 65:- y £ + 6 & ( 6) 

etc. Note that the Eulerian bias bk depends on the Lagrangian bias factors of equal and lower order, and the overall 
structure is precisely that shown as the monopole contribution to the bias in Q . 

Accounting for the fact that the evolution of the environment will, in general, be nonlocal means that one simply 
replaces the 8q(8) mapping with the nonlocal one. It is straightforward to check that carrying this through, order by 



3 



order, yields the additional nonlocal bias terms given in [4. Il5|. However, if halo bias is nonlocal even in Lagrangian 
space, then this will provide additional contributions to the Eulerian bias factors. To see the structure of these terms, 
it is useful to consider models of halo formation in which factors other than the initial overdensity are important. 



B. The triaxial collapse model 

In triaxial collapse models, e.g. [l2j], the evolution of a patch is determined by more than just its internal overdensity. 
The simplest of these models uses the fact that, at each position in the initial field, one may define the deformation 
tensor, D, whose elements are the second derivatives of the gravitational potential. The overdensity, which is the only 
quantity which matters for the spherical model, is the trace of this 3x3 matrix. So the question arises as to which 
(combinations) of the other elements of this matrix matter? 

To describe the shape of the gravitational potential it is common to introduce the ellipticity e and prolateness p, 
defined from the eigenvalues Aj (i = 1, 2, 3) of Vy 

\ , \ , \ Ai — A 3 Ai + A 3 — 2A 2 , , 

d = Ai+A 2 + A 3 , e= — — — , and p= — . (7) 

This set of parameters is used in triaxial evolution models of nonlinear structure formation [l2l.[l3l]. However, because 
e and p are ratios of the eigenvalues, it is not obvious that they are the best choice of parameters in a perturbative 
analysis. In particular, one might have wondered if the rotationally invariariant quantities, 

h = Tr(£)) = ^;A i = «5o, (8) 

i 

I 2 = AjA 2 + AiA 3 + A 2 A 3 , (9) 
h = Det(D) = JjAj (10) 

i 

are more relevant. When expressed in terms of (5q, e,p) these are 

h = I [1 - (3e 2 + P % h = |(1 - 2p)[(l +pf 9e 2 ] (11) 

Since the Ij do not depend on taking ratios of the eigenvalues, they, or other quantities built from them, have 
considerable appeal. Indeed, G2 = — 2/ 2 and £3 = 6/3 are the fundamental quantities in Q. 
Another interesting combination is 

2 r 2 or ,2 fo 2 , 2, , , _ 21? - 9hl 2 + 27/3 _ 2S*p(9e 2 -p 2 ) 



S Q = I lt q£ = If - 3I 2 = 5$ (3e 4 + p*) , and 
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Despite the appearance of I\ = Sq in their definition, q and u are actually independent of 5$. Moreover, they are 
precisely the quantities which arise in a perturbative analysis of the ellipsoidal collapse model: 3\ and J 2 of [16| are 
our 5q and 9mq respectively. A final combination which also arises in triaxial collapse models is 

v = Ai — A 2 and w = S — 3A3, (12) 

where < v < w and < w < 00 [l7j . Like q 2 and u 3 , v and w are also independent of 8. 

Notice that e, p, v, w, q 2 , and Uq all vanish for a spherically symmetric perturbation (Ai = A 2 = A3 = 5q/3). But 
in general, qo follows a xi distribution [l8| (they actually called it r; we have changed notation to emphasize the fact 
that q 2 is the quadrupole which arises in perturbation theory analyses). 



C. Nonlocality from evolution in such models 



Whereas the spherical evolution model assumes that halos are associated with regions where Sq > S c , triaxial 
evolution models generically assume that Sq > S c (e,p) or S c (qo,uo); the critical density required for collapse depends 
on the other elements of the deformation tensor. Most work to date has treated the effect of (e,p) approximately, by 
using appropriately defined average values to estimate halo abundances [l3j and how they correlate with the large scale 
environment [l8|. In this approximation, Lagrangian halo bias remains local, but the bias coefficients are modified 
compared to the spherical case. In addition, halo bias at late times is also treated approximately, by assuming that 
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the large scale environment evolves according to the spherical model. As a result, Eulcrian halo bias is also local 
(see [H[[i(| for explicit expressions). If one keeps the local approximation for Lagrangian bias, but accounts for the 
nonlocality of nonlinear evolution of the dark matter, then this will again yield the same structure for the nonlocal 
terms that was described in and shown to be robust beyond the assumption of conserved tracers. 

Before reconsidering the question of local Lagrangian bias, it is worth noting that the nonlocality of nonlinear 
evolution is part and parcel of the ellipsoidal collapse model. This is because, in this model, the second and third 
order approximations to the nonlinear density are given by [l6j 

*» = + (13) 

, 3) 341 r3 338 . 2 92 3 

S { ' = <5n H <W H u?,. (14) 

567 945 m 441 v ' 

Comparison with Eqs. (34)- (36) of [|| shows that the terms proportional to 5q and gjj in the expression for 5^ 

(2) (2) I — 73t 

are essentially the same as their monopole fC\ /_ and quadrupolc K\ /_ 2 (also see [22|). Differences between this 

approximation and the exact dynamics appear for where the monopole is the same as lC^j =Q , but the quadrupole 
differs. In addition, there arc differences coming from nonlocal potentials that show up at third order (see e.g. 
Eq. 107 in 0), and dipole terms that correct both second and third order expressions. This shows explicitly that 
a self-consistent use of the triaxial evolution model should yield a better description of halo bias. In particular, if 
one uses this model for the mapping between Lagrangian and Eulcrian bias, then it makes sense to reconsider the 
assumption that the bias is local in Lagrangian space. 



D. Nonlocal Lagrangian bias for triaxial collapse 

Triaxial evolution models will give rise to nonlocal Lagrangian bias if the parameters which determine S c , e.g. (q, u) 
or (e,p), couple to the large scale environment. To see that this is generic, consider a simple nonlocal model in which 
halo abundances depend on the tracelcss part of the initial shear field qo as well as on the initial local density So, an d 
that this arises because the critical density required for collapse depends on the shear field. This model is particularly 
simple because (So/cro) an d (<7o/co) 2 are independent, the former being drawn from a Gaussian variate with unit 
variance and the latter from a chi-squared distribution with 5 degrees of freedom (l8| . This simplifies the analysis 
required to estimate the numerator and denominator of 

1 + S h (m\6 , q ) = = 1 + b 1 d Q + b 2 — + c 2 + . . . (15) 

n(m)Vo 2 2 

In general, notice that if we integrate over all values of q$ at a given Sq, then this will yield 

l + 8 h {m\5a) = l + %6 + $$ + <%y^- + ... (16) 

Since (<Zo|<5o) = (<7o) does not depend on So, non-zero values of c 2 will make it appear as though the bo component 
differs from unity. And since (q^) = ctq, this offset from unity will be vanishingly small on large scales. 

E. Nonlocal Eulerian bias for triaxial collapse models 

Combining this nonlocal Lagrangian bias with Eq. (| 14[) for nonlocal gravitational evolution yields 

1 + = (l + S)(l + %5 + %$ + ( ^ q l + ..^=i + b\8 Q + b^ + ^ + 5 + b\5oS 

= i + 5{b\ + l)+^b\/2l + b\) + ^{c\-Sb\/2l). (17) 
Since, to lowest order, q$ = q 2 , this makes the Eulerian bias factors 

h = l + %, b 2 = b\ + ^-b\ and C2=C L-A^ (18) 
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Comparison with equation (|6]) shows that b\ and b 2 are related to the Lagrangian bias factors just as in the local 
spherical collapse model; the nonlocality shows up as a nonzero value of c 2 . 

To make the connection to 3 we now express this in terms of Q 2 = (2/3) (qq — <5q). This makes 

1 + 5 h (m\8o, &) = 1 + b\ So + Q$ + (4/3) 72 L ) f + 7 2 L Q2 + ■ ■ ■ , (19) 

where we have defined 72 = 3c 2 /4. Since Q 2 = (2/3)(gg-d§), we have that (g 2 \S ) = (2/3)((gg)-<Jg) = (2/3)(c^-5§) -> 
— (2/3) Sq on large scales. So, 

1 + S h (m\S ) = l + b\So + b\^ + (4 72 L /3)^ .... (20) 



Similarly, 



2 v / 2 



1 + = 1 + ^(^ + 1) + y (8#/21 + #) + ^±|/^( c l -86*721) 

= 1 + + 1) + y (6& + 4 72 L /3) + S 2 ( 72 L - 2#/7), (21) 



making the Eulerian bias factors 



h = l + b\, b 2 = ^b\ + b\ + ^- and 72=7 2 -^i- (22) 

If averaging the distribution of Q 2 at fixed S yields the same as in Lagrangian space, i.e., (Q 2 \S) — — 25 2 /3 (and we 
set (To — > 0) then equation ([2Tj) implies that 

1 + 5 h {m\5) = 1 + b, 5 + (b 2 - 4 72 /3) 6 - + . . . (23) 

which is consistent with Eq. (117) of And, when c 2 = 0, then our 72 = —25^/7, which is consistent with Eq. (118) 
of 0. 

In the next section we use the excursion set approach to estimate the numerator and denominator of Eq. (|15p . i.e. 
(N m \5o, <7q) and n(m). Expanding in a Taylor series yields predictions for the bias factors b\ and (or 7^), and 
hence for the Eulerian bias factors. 



F. Cross-correlations 



The expressions above imply that the Lagrangian space cross-correlation between halos and mass is 

(* r [1 + 5 h (m\5 , ql)}) « b\ (S r S ) + b f(S r S 2 Q ) + ^(5 r q 2 ) + |(<5 r <5 3 ) + . . . (24) 

where the average is over the joint distribution of So and qo at one position, and of So at another position a distance r 
away. This means that the term (5 r qo) should be thought of as (S r (qQ\So)) , where the inner average is over values of 
qo at fixed Sq, and the other average is over all S r So pairs separated by r. This shows that replacing Sf,,(m\So, qfi) by 
<5^(m|i5o), its mean value for given So, before measuring the cross-correlation should yield the same answer as if one 
had included the full scatter. This explains the agreement between no-scatter and full-scatter measurements presented 



in Figure 1 of [15 1. 



For Gaussian initial conditions, all terms of the form (S r 5o) in the expression above can be written as (5o(S r \So}} = 
(<5q +1 ) (S r So) / (Sq) ■ Thus, the entire local bias contribution is linearly proportional to £(r) = (S r So) [ll|. Of course, in 
the present example, q$ and So are independent, so all terms of the form (<5 r ?o ) vanish, meaning there is no nonlocal 
contribution to the cross correlation. On the other hand, the auto correlation will receive contributions from terms 
of form (S^qo k ). These will first appear at order (S^qo) and (qfqo); since they are of the same order as (S^Sq), they 
will also contribute to the bispectrum. 
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III. NONLOCAL LAGRANGIAN BIAS IN THE EXCURSION SET APPROACH 



The analysis above is useful but otherwise empty formalism. The main goal of this section is to see if the nonlocal 
bias factors are comparable in magnitude to the usual ones, by estimating how they depend on halo mass and time. 
We use the excursion set approach to do this. 

In the excursion set approach, the abundance of objects of mass m is closely related to the distribution of scales R m 
on which the initial overdensity S > B for the first time, where B ^ S c is a measure of how the 'barrier' for collapse 
differs from the value 5 C for the spherical evolution model. In triaxial collapse models, B is a function of e and p 
[lH, but, in what follows, we will use a simpler model in which B is a function of q 2 (Eg. IT21 shows q 2 is a particular 
combination of e and p). We emphasize that the logic which underlies our approach is not confined to this choice; 
we have chosen this simple model because it provides a particularly easy way to see the main effects associated with 
non-spherical collapse. 

The particular model we will study has 

B(q)=8 c {l + y/#M), (25) 

where q c is a characteristic scale for the effects associated with the shear. This expression is motivated by setting 
p = in Eq. (3) of (l3j |: with q 2 w 85 2 it provides a reasonably good (but by no means perfect) description of the 
actual initial overdensities in patches which form halos. Note that non-zero q always yields B > 5 C , a property which 
will play an important role in what follows. 

The dependence of B on q means that we must construct a list of 5i and qt which exhibit the correct correlations 
between steps, and then find the smallest m (largest smoothing scale) for which 8 m > B(q m ). (We have deliberately 
used m as the index for the smoothing scale, since this scale is monotonically related to the inverse of the mass.) 
This is straightforward to address using Monte Carlo methods @, HH, but, as we describe below, it turns out to be 
possible to write down rather accurate analytic approximations to the main results. We will first consider the case 
in which the steps in the walks are uncorrelated, before describing the modifications which come from accounting for 
correlations. 



A. Universality associated with uncorrelated steps 



While the Appendix describes a more careful analytic calculation for walks with correlated steps, we have found 
that the substantially simpler analysis outlined by |l8| actually provides a rather good approximation to the results 
for walks with uncorrelated steps, and yields considerable insight. This analysis replaces q 2 with its mean value a 2 . 
This means that B(q) — > B(a), and so the first crossing distribution associated with the six-dimensional walk (one 
for 5 and 5 for q 2 ) should be well approximated by that for the one-dimensional condition S > B(a 2 ). Reasonably 
accurate models for the first crossing distribution of such moving barriers are given in (l8j . In particular, if 

%) = Mi + \^M)^Mi + \AW), (26) 

then the first crossing distribution is well-approximated by 



/ 1/4 \ P -("+\/^I) 2 /2 

vj{v) a 2 [v + -J= - = where v = 8 c /a. (27) 

\ \A?c/3? / V 2tt 

In this model, the ratio S c /q c is a measure of the strength of the nonlocal effects; these disappear in the limit S c /q c — > 0. 

For similar reasons, the first crossing distribution for walks which start from some non-zero 6 = So and q 2 = q§ 
should be well-approximated by that for one-dimensional walks which must cross 



B(q)^S c [l + ^(a 2 +q 2 )/q 2 ] - S . (28) 

The term in the square root follows from the same logic as for the unconditioned walks; i.e., one replaces the 
dependence on q 2 by a dependence on the mean value (q 2 \qo), where the constraint go means that q is now drawn 
from a non-central chi-squared distribution (with 5 degrees of freedom). This shows that the constraint enters with 
a plus rather than a minus sign. 

Therefore, the local bias parameters are given by the usual derivatives with respect to S c , whereas the nonlocal bias 
parameter is slightly more involved. In particular, the result of using this barrier when estimating the first crossing 
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FIG. 1: Excursion set description of the nonlocal Lagrangian bias which results from the local shear affecting the collapse 
threshold (Eq.[25]with qj: — 85%). Left: Distribution of first crossing scales for unconditioned walks (black); walks which began 
from 5q/8 c = 0.2 but qo = (red); walks which began from (qo/qc) 2 = 0.4 but So — (magenta). Dotted curves show the 
Press-Schechter and Sheth-Tormen distributions, and solid curve shows our Eq. (|27[) which follows from approximating the 
collapse barrier following Eq. (126p . Right: Associated large scale bias factors; dotted curves show our predictions (Eqs. 1291 
and 130)1 . 



distribution, expanding to lowest order in So and q^, dividing by the distribution associated with (5q, qo) = (0, 0) and 
subtracting one, yields: 



6 C b\ = v 2 - 1 



1 



5 2 c b\ = vH z {v) 



q c /S c l + Auq c /8 c 



2v 2 



q c /S c l + 4vq c /5 c 



- <? r L 



and 



X 2 L 

0„ Cn = 



(qc/5cY 



1 + v(q c /S c ) - 



%q c /5 c f 



1 + (128/35)m? c /<5 c 



(29) 



(30) 



Figure [T] shows that this works rather well. The plot shows results for q 2 /S 2 =8 for reasons that will be clear shortly, 
but we have checked that it works well for a wide range of q c /S c - The panel on the left shows the distribution of first 
crossing scales estimated from 64000 walks: black, red and magenta are for walks which start from (So/S c ,q 2 /5 2 ) = 
(0,0), (0.2,0) and (0,0.4), respectively. The first represents the unconditional distribution, whereas the second and 
third were chosen to isolate the effect of a large scale overdensity and (quadrupole) shear, respectively. Notice that 
the model predicts more massive halos in regions with positive So, but fewer in regions with non-zero qo. This is a 
direct consequence of assuming that the barrier height is an increasing rather than decreasing function of q 2 . I.e., 
when v S> 1 then equation (|2"B")) becomes B/a — > v (1 + qo/q c — Sq/5 c ). In the opposite limit, v <C (5 C / q c ) / (qo/ q c ), 
B/a — > v + (S c /q c ) — <5 /cr; for sufficiently low mass halos, only Sq matters. 

The solid curve shows Eq. (|27|) ; it describes the unconditional distribution (i.e. black symbols) well. The dotted 
curves show the Press-Schechter and Sheth-Tormen distributions; for ST, we have set a = 1 (rather than 0.7). Since 
the shape of the ST curve provides a good description of real halos, the good match between the solid and ST indicates 
that, for this choice of q c /S c , the shape of the halo mass function is like that in data. In turn, this suggests that our 
model should yield realistic estimates of the sign and amplitude of nonlocal effects. For example, this match indicates 
that the trend to have more massive halos in regions with positive 5q, but fewer in regions with non-zero q 2 is realistic. 

The panel on the right shows the ratio of the environment-dependent distributions to the unconditional one, scaled 
by (<W^c) an d {qo/qc) 2 respectively. The dotted curves show 5 c b\ and 5 2 of Eqs. (|2T))) and (|3"U)) ; they describe 
the measurements rather well, suggesting that our analytic treatment, which accounts for the stochasticity in the 
barrier distribution by ignoring it following [l8| , has captured the essence of the problem. We show elsewhere how to 
modify this trick for dealing with stochasticity so that it also works for other barrier shapes. Note in particular that 
our analysis indicates the magnitude of can be comparable to that of b\ , so the Eulerian nonlocal bias parameter 



C2 



86^/21 can be substantial, particularly when v 1 (the most massive halos). 



8 




o 1 
ln(y = 6l/a s ) 



FIG. 2: Comparison of excursion set predictions for the first crossing distribution (left) and associated nonlocal Lagrangian 
bias factors (right) for a flat ACDM cosmological model, when correlations between steps have been included in the analysis. 
Panel on the left shows results for walks with uncorrelated steps (labeled sharp-k), and for walks in which correlations arise 
from TopHat smoothing. The corresponding first crossing distributions are well-described by Eqs. (|27[) and (|31[) . respectively 
(solid lines). (The curve showing Eq. (|A4|) is almost indistinguishable from that for (|31|l . so we have not bothered to show it.) 
We have only shown the bias factors for the TopHat smoothing filter; they are well-described by the solid lines which show 
Eqs. (p3)l and (p4)) . 



B. Departures from universality from correlations between steps 



The analysis above was based on walks with uncorrelated steps. These correspond to smoothing the initial Gaussian 
field with a filter which is sharp in fc-space. Smoothing with filters which are more localized in real space (e.g. a 
spherical tophat) will result in walks with correlated steps; because such smoothing filters are intuitively closer to the 
physics of collapse, predictions which are based on correlated steps are expected to be more realistic (see [IJ HE HE] 
for why this is not the full story). 

In what follows, we will show that accounting for such correlations turns out to be relatively simple, and has some 
important consequences. The analysis is simplified for two reasons: S and q 2 are independent whatever the smoothing 
filter, and the argument about transforming the six-dimensional walk problem to an effective one-dimensional barrier 
should work even if steps are correlated. Recently, accurate models for the first crossing of moving barriers by walks 
with correlated steps have become available [E] (their Eq.5), so, combining these for the square-root barrier problem 
that is relevant here, yields 



/2tt 



erfc(IV/\/2) , e - r " / 2 

1 — -l 



2-kTv 



where v = ^ and r 2 = J**^ , (31) 



where S' = dS/da 2 . For ACDM, T 2 w 1/3. 

The Appendix contrasts this with the approximation suggested in Q , in which the first crossing distribution should 
be thought of as averaging the distribution for fixed q over the distribution of q. At large v 3> 1 both approximations 
predict a factor of 2 fewer objects than the solution for uncorrelated steps; this difference decreases as v decreases 
until sufficiently small v <C (l/4)/(g c /£ c ), when the T-dependent factor begins to dominate. Since v ~ 0.1 is below 
the regime of most cosmological interest, and it is also in the limit where the analytic approximation of Q is expected 
to break down anyway, we will restrict attention to larger v. In this regime the first crossing distribution will differ 
only slightly from that for uncorrelated steps, so we expect to find similar bias factors with the following caveats. 

First, since T depends on the shape of the power spectrum, we expect plots of, e.g. C2 versus b\ to no longer be 
universal, but to depend on P(k). However, because T is defined by a ratio, changes to the overall normalization of 
P(k) will cancel out. Second, the bias factors b n become fc-dependent because of the correlation with the curvature 
term 5' If q were correlated with 5' then C2 would also become fc-dependent. However, q is independent of 6' , and, 
in any case, in what follows we will restrict attention to the scales on which the fc-dependence can be ignored. 
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FIG. 3: Dependence of excursion set prediction for the relation between Si and b\ on the shape of the power-spectrum. The 
two solid lines are for correlated steps with T 2 = 1/3 (lower) and 2/3 (upper); dotted line is for uncorrelated steps and is 
independent of P(k). 



For walks with correlated steps that are constrained to pass through (Sq, qo) on scale So, the first crossing distribution 
can be derived similarly to when the walks started from (0, 0) on scale So = 0. The analog of equation (|28p becomes 



B(q) 



1 + J [a + (S x /S ) 2 (q 2 - So)]/q 2 - (S x /S )<5 , 



where 



Sx 
So 



(SSp) 
(S 2 o) 



(32) 



When v ^> 1 and So — >• then B v[l + (S x / So){qo/qc — So/S c )]' except for the factor of S x /So, this is the same 
shift as for walks with uncorrelated steps (compare Eq. I28|) , so we generically expect the same qualitative trends with 
environment. 

Expanding the analog of equation (|31j) in powers of 5o and qo-, and taking the So limit yields 



6 C b\ = v 2 - 1 



q c A 



+ A(Tu), 



vH z {v) 



and 



r2 L 

°c c 2 = - 



(qc/5 c y L 



1 + v{q c /5 c ) 



| 2isH 2 (is) 

q c /S c 

2(q c /S c ) 



(q c /s c y 



+ (T 2 v 2 + 2) A(Tv) 



2v 

q c /5i 



A(Tu), 



where we have defined 



A(x) 



1 + crf(x/V2) 



-x 2 /2 



(33) 



(34) 



(35) 



The Appendix describes the large scale bias associated with equation (|A4|) . 

Figure [5] compares these predictions with Monte-Carlos, showing that the model works rather well. Following (To| . 
we estimate the bias factors in the Monte-Carlos using cross-correlations. Namely, we estimate 



h = 



[sf(s)} \fSo~ (Sx/So) 



b 2 = 



( S f( S \So)(S 2 /S 1)) 
( S f( S ))S (S x /So) 2 ' 



C2 = 



(sf(s\So)(r 2 /(r 2 )-l)} 
[sf( S ))V a r(r 2 /(r 2 ))(r 2 )(S x /So) 2 ' 



(36) 



where the angle brackets denote sums over walks which first cross at s, V&i^q/ (r 2 )) denotes the variance of r 2 1 (r 2 ) 
on scale So, and the factors of S x /So must be included when making estimates in this way for reasons given in 1C|. 
We present results for S = 0.09, for which S x /S ~ 1.45. 

The analysis above shows that, as was the case for uncorrelated steps, b±, b 2 and c 2 are all expected to depend 
on v. However, there is now an additional dependence on the parameter T. As a result, for walks with uncorrelated 
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steps, there is a predicted relationship between C2 and b\ which is expected to be independent of the underlying 
cosmological parameters or power spectrum. Correlations between steps introduce departures from this universality, 
so it is interesting to quantify this dependence. Figure |3j shows that the predicted departures are small. Thus, 
our analysis indicates that our main finding - that the Eulerian nonlocal bias parameter C2 = — 86^/21 can be 
substantial for the most massive halos - is robust to changes in the power spectrum. 

IV. DISCUSSION 

If halo bias is nonlocal in Lagrangian space fEq. [T5|) . then this will add nonlocality in Eulerian space bias (Eg. H8|) . 
We provided an explicit calculation of this effect, in which halo formation depends on the initial local density So and 
shear field q^. 

In the excursion set approach, the problem of estimating halo abundances reduces to solving for the first cross- 
ing distribution of a suitably chosen barrier (Eq. I25|) by 6-dimensional random walks. We provided a simple but 
accurate analytic approximation for this distribution for the case in which walks have uncorrclated steps (Eq. 1271) ; 
the approximation follows from treating the full 6-dimcnsional problem as an effective one-dimensional one (Eq. 1261) . 
Predictions for halo bias come from studying walks which do not start from the origin. We argued that the associated 
first crossing distribution is best thought of in terms of a shifted barrier (Eq. |2"5| , from which it is straightforward to 
derive formulae for halo bias formulae (Eqs. |2T)1 and I5D"| . These formulae, which quantify how the large scale density 
and shear fields affect halo abundances as a function of mass and time, are quite accurate (Fig. [1}. 

For walks with correlated steps, the predicted first crossing distributions and bias formulae can be written in units 
in which they are universal (independent of power spectrum and cosmology) . We argued that this universality should 
be weakly broken if steps are correlated (Fig. [3]). In this case, too, we provided analytic approximations for the 
unconditional first crossing distribution fEq. I3T|) and halo bias factors fEqs. 1551 and IM|). which were quite accurate 
(Fig. [2]). Our results indicate that nonlocal bias effects, as quantified by the parameter C2 = — 86^/21 of Eq. (|18l) . 
can be substantial for the most massive halos. 

Our analysis is easily extended to describe the more complicated case in which the barrier S c de pen ds on (e,p) of 
Eq. (JT]) rather than simply q 2 . This is because £ e (e, p) is a function of the combination eS and p5 [13, and, like q 2 , 
these combinations are actually independent of 5 [l8j ] - Alternative parametrizations of this (l7T ] have S c (v 7 w) where v 
and w (defined in Eq. fT2|) are also independent of 5. This means that the analog of Eq. (|A4|) . in which one averages 
over a distribution of first crossing distributions, remains a good approximation. 

In addition, the idea that one can map the 6-dimensional walk problem to a 1-dimensional moving barrier problem 
should continue to hold when 5 c {e,p) or <5 c (v, w) rather than 6 c (q), so the analog of Eq. (|31|) should also provide 
a reasonable approximation. Therefore, we believe our analysis should be applicable to other parametrizations or 
models of the effects of nonlocality. Indeed, our analysis should also apply to cases where one places conditions on 
the eigenvalues of the deformation tensor (e.g. all three have the same sign), rather than the combinations e,p, or 
v,w. In this respect, it provides the basis for modelling not just halos, but the abundance and spatial distribution 
of supcrclusters, filaments, sheets and voids as well. In particular, our analysis predicts that all of these constituents 
of the cosmic web should exhibit nonlocal bias effects; it will be interesting to see if such effects are discovered in 
simulations. 

In our analysis of nonlocal halo bias, we assumed that the effect of the large scale environment was to affect the 
distribution of 6 and q on smaller scales, but not the s hap e of the collapse barrier. If the halo formation process depends 
on the surrounding environment, as the analysis in |25[ suggests, then this will provide an additional contribution to 
nonlocal bias. It is straightforward to include such an effect in our treatment of conditional n-dimensional walks, but 
our current expressions do not do so. 

We stated at the start that the nonlocal effects which are the subject of this paper, and which enter at the quadratic 
level (q 2 , e<5, etc.), are qualitatively different from those which enter even at the linear level, and contribute fc-dependent 
terms to the bias. In principle, both effects are present in our correlated walks calculation - the latter arise from the 
dependence of the first crossing distribution on the derivative of the initial density field [1, [l(| , rather than on the 
anisotropic distribution of the mass. Since our main goal was to illustrate the effects of q 2 , etc., we ignored these 
other terms, but a more complete model would include them. Separating these terms from one another should be 
possible, since the terms which depend on derivatives of the field will only contribute to the monopole of the bias. 
This, and comparing to measurements in numerical simulations, is the subject of work in progress. 

Quantitative agreement with simulations may require merging the formalism here with the additional requirement 
that one is interested in special positions (such as peaks), rather than random positions in the field. As noted in 
[Tlj . the formalism of [|| on which our analysis is based allows the inclusion of this requirement with no additional 
conceptual complications; for peaks in a Gaussian field, this is particularly straightforward, because the distribution of 
the shear field around such peaks is known [26[ . Our analysis also suggests that a fruitful extension of the peaks model 
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to smaller masses than that on which it usually breaks down is to look for peaks in the field defined by 5 — \/ q 2 jq 2 . 
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Appendix A: Analytic estimate of nonlocal Lagrangian bias 

For barriers of the form 5 C + q5 c /q c , Eq. (13) of Q yields a simple estimate of the first crossing distribution: 

fc(Ti//V2) ver^^l 2 



yfiy 

where 



2-kTv 



e -(^+Vp 2 ^/? c 2 ) 2 /2 

dpp 5 (p) -j= , (Al) 



P 2 \ 2 ) r(5/2) 



p 5 (p) dp = — [ — ) -^7777^, with p=q/<7- (A2) 



The integral can be written in terms of the parabolic cylinder function U(a, z) with a = 9/2 and z = (1 + 5q 2 / S 2 ) x l 2 , 
which, in turn, can be written in terms of derivatives of e z / 2 erfc(z/-\/2)- The integral equals 

.e- v2 ' 2 ( 5/2 \ 5/2 4! //4rr/n/0 .,^ V/2 r(l/2)/ 5/2 \ 5/2 d*[e*'/*eTfc(z/y/2)} 



'^U/^ + sJ r(5/2) gZ ^ (9/2 ' z) - ^ F(5/2)Uc 2 /9c 2 + 5j dz± (A3) 
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FIG. 4: Dependence of excursion set prediction for the first crossing distribution of the barrier S c (1 + q/q c ) by six-dimensional 
walks, on q c . Solid curves show Eq. (1A4[) . and dashed curves show the approximation Eq. (|31[) . Curve labelled MS shows the 
q c — > oo limit, and dotted curve labelled BCEK shows this limit for one-dimensional walks with uncorrelated steps. For the 
walks with correlated steps, the underlying power-spectrum is ACDM, and the correlations are due to smoothing with a real 
space Tophat. 



where z = (1 + 5^/^)-i/2 ThuSj 



"/to 

c 2 



D -,s 2 /2 



1 



erfc(r^/\/2) 



V2tt 

(5 + (<5 c / 9c ) 2 ) 4 
2 (u5 c /q c ) 
'(5 + (<W<7c) 2 ) 3 "' (5 + (S c /q c W 



2 V2Trzy 

6 (v5Jq c ) 2 



(5/2) 5 / 2 
r(5/2) 

3 



[Ci(i/) - 



where 



2tt 



(5 + (<5 c / 9c ) 2 ) 3 (5 + (<5 c /g c ) 2 ) 2 
2(vS c /q c ) 3 



5 + ((5 c /g c ) ; 



r erfc(j/)e 2 ' 2/2 



(A4) 



Although this expression is not particularly illuminating, we have included it because the peak background split bias 
factors are easily estimated from the fact that 



dm[eZ2 'Z {a ' Z)] = (-I)" (V2 + a) m U(a + m, Z ). 
dz m 



(A5) 



More insight comes from approximating the integral over p in Eq. (|A4[) with the value of the integrand at its mean 
value (p) = 1. This yields Eq. (f31"T) in the main text. A little algebra shows that this particularly simple and intuitive 
approximation should be quite accurate when q c /S c 3> 1. Figure 0] shows that this is indeed the case, although 
equation (|A4[) is accurate over a wider range of q c /S c . 

To estimate the large scale bias, for comparison with the peak background split estimate, we must first estimate 
the first crossing distribution for walks which pass through (A, Q) on some S < s. Then the Lagrangian bias factors 
are given by expanding 



l + S h {m\A,Q) 



= f{m\A, Q) 
f(m) 



(A6) 



in powers of A and Q. This shows that one generally expects halo abundances and hence Lagrangian bias to depend 
on Q as well as A. 

The effect of the constraint is to modify the distribution of (5, q) but not the barrier (this is actually an assumption 
- in principle, the shape of the barrier could depend on the large scale environment, in which case this would provide 
an additional environmental effect), so the calculation is actually rather similar to that for unconstrained walks. The 



13 



independence of q and S on all scales means that the conditional distribution, given the values A, Q on some other 
scale (in what follows, we will assume this other scale is larger), factorizes: 



where 



ff(*|A) 



cxp 



piS', 6, q\A, Q) = P5 (q\Q) g(S'\6, A) g(6\A), 

.,2 <^) 2 



-(5- 7V ' S /SA) 2 /2 S (l-7 2 ) 



^27^(1-7*) 
and ps(r\R) is a non-central chi-squared distribution: 



with 



(<5 2 )(A 2 



(SAy 
sS 



(A7) 



(A8) 



Ps(q\Q)dq= d Pl p 5 ( Pl )e "d-V) 



5 7 2 (Q 2 /S) 5p 2 \ T(5/2) 
2(1 - 7 2 ) 2 J r(5/2 + j) 



where 



g> 2 
1 - 7 2 



(A9) 



Notice that this is simply a %| distribution for Pl when Q = 0. In general, the integral over p 7 can again be written 
in terms of parabolic cylinder functions. But, to lowest order in Q, 



Pb{q\Q)dq = dp 1 p b {p 1 ) l-(l-p ) 



, 5(7 2 Q 2 /g) 
2(1 -7 2 ) 



where 



2 _ g> 2 

Pr > ~ 1 - 7 2 



(A10) 



This is most easily understood by thinking of the non-central chi-squared distribution as a Poisson mixture of central 
chi-squared distributions of ever higher order, and keeping only the lowest two terms. Note that, for a sharp fc-filtcr, 
7 2 = (A 2 )/((5 2 ) = S/s, making j 2 R 2 /S = R 2 /s in the expression above. However, we are interested in the general 
case, for which it is convenient to define S x = (6 A) making j 2 R 2 /S = (S* x / S) 2 (R 2 /s). 



Vb{q\Q)dq = d Pj p 5 (pj)e Md-y) 



d Pl p b { Pl )e Md-V) 



E 

3=0 

E 



5 7 2 (Q 2 /s)Y T(5/2) 



E 



2(1 - 7 2 ) J r(5/2+j) 

5 7 2 (g 2 /5) a; v i ^/57 2 (g 2 /^ 



2(1 -7 2 ) 



E 

3>i 



i=0 

lira /c\..\ j-» 



2(1 -7 2 ) 



£(5/2) 
r(5/2 + i)VC?-i)!i! 



If we define 



6-jy/I/SA (5- 7 y77^A 
V2s(l-7 2 ) ^ 



on large scales. Then 



erfc(r^ A /V2) e~ r2 ^/ 2 



+ 



27T IV A 



dqPb(q\Q) 



.-(^A+vVif/gf) 2 ^ 



(All) 



(A12) 



(A13) 



The purely local contributions to the bias come from considering the ratio f{v\A,Q = Q)/f{y) in the S — > limit. 
These correspond to replacing ^ — > in equation (|A4|) . In addition, the expression above shows that if the barrier 
does not depend on q, then vf{v\A,Q) does not depend on Q; nonlocal effects on the bias are entirely due to 
dependence of the collapse barrier on q. 



